Spherical non-linear absorption of cosmological scalar fields onto a black hole 
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In this paper we track the non-linear spherical evolution of a massless scalar field onto a 
Schwarzschild black hole space-time as a first approximation to the accretion of cosmologically 
motivated classical scalar fields. We perform an analysis related to wave packets described by wave 
number and width. We study various values of the wave number and found that for k = and width 
packets bigger than the Schwarzschild radius, the absorption is not total. In the cases we studied for 
k > 0, the black hole absorbs the total amount of energy density of the scalar field moving toward 
the horizon. Our results indicate that assuming spherical symmetry, in the non-linear regime, there 
are cases for which scalar fields are allowed to survive outside black holes and may eventually have 
life-times consistent with cosmological time scales. 

PACS numbers: 04.40.-b,04.25.D-,95.35.+d,95.36.+x 
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I. INTRODUCTION 

Cosmological scalar fields have been an important in- 
gredient of cosmological models since the model of infla- 
tion proposed a mechanism providing exponential cosmo- 
logical growth [l|. Considering the supernovae redshift 
observations a smooth classical scalar field was proposed 
to play the role of dark energy Q . And scalar fields have 
also been proposed to play the role of dark matter, first 
at galactic scale [|| and then at cosmic scale On the 
other hand, current observations indicate that there are 
supermassive black holes in the centers of galaxies with 
masses of the order of 10 6 — 10 9 solar masses One of 
the quests consists in determining the components feed- 
ing the black hole, mainly related to the accretion of ei- 
ther baryons or dark matter. For instance it can be con- 
sidered that supermassive black holes grow through the 
accretion of baryonic and dark matter onto seed black 
holes of intermediate mass of about 10 3 to 10 4 M Q Q. 
Also some results indicate that only about 10 % of their 
mass is due to the accretion of dark matter (e.g. @), and 
more mature models involving the study of the Newto- 
nian phase space [8| indicate that the time-scale for the 
accretion of collisionless dark matter to contribute signif- 
icantly to a supermassive black hole mass is too long (see 
e.g. |9j]). Currently the study of ideal gas dark matter 
accretion onto supermassive black holes is under research 
and important bounds on the equation of state of dark 
matter may arise [To|, EH ■ 

Thus studying the accretion of cosmologically moti- 
vated scalar fields is expected to have an astrophysical 
impact within the subject of the viability of dark mat- 
ter candidates. In fact the accretion of scalar field dark 
matter and dark energy onto supermassive black holes 
has been already explored under certain assumptions. In 
[ijj the accretion rate of scalar fields is calculated in 
order to estimate whether or not the accretion rate is 
consistent with the mass and life time scale of supermas- 
sive black holes. Later on, in [l3| the same problem was 
treated considering that the scalar field has two ends, the 



black hole's horizon and future null infinity, which was 
achieved by using hyperboloidal slices of space-time and 
compact radial coordinate [l4| which allowed the scalar 
field not only to be accreted by the black hole, but also 
to leak through scri+, and it was found an extremely 
high dilution rate of the scalar field density that may 
rule our scalar fields as dark matter and dark energy On 
the other hand, more recently in [l5| it was studied the 
accretion rate in terms of the wave number of a wave 
packet of scalar field, and found interesting results re- 
lated to dependence of accretion rates on the width and 
wave number of wave packets, that is, for wave packets 
with small wave number and very spread the absorption 
of the scalar field is lower than for thin packets with big 
wave numbers; this problem was also studied under a dif- 
ferent context related to the absorption of light in 
In 17] the same problem was treated and it was found 



that there are particular initial scalar field profiles that 
allow long-lived scalar field densities around a black hole. 
All these results have been obtained considering spheri- 
cal symmetry and that the scalar field is a test field, that 
is, the geometry of the space-time remains fixed. 

In this paper explore a step forward and include the 
evolution of the space-time geometry. That is, we study 
a spherically symmetric scalar field evolving onto an ex- 
isting black hole space-time, considering the evolution 
of the geometry too, so that we can measure the black 
hole's horizon growth and determine more precise accre- 
tion rates and absorption ratios. We focus on the study 
of the absorption of the scalar field in terms of the pa- 
rameters of wave number and width of the initial wave 
packet as described in (l5j . 

In order to clearly set the astrophysical scenario we 
deal with, we consider important to stress that the black 
hole is already formed initially, and that it is asymptoti- 
cally flat. These conditions have important implications: 
(i) we assume the cosmological effects on the system are 
neglected, and the evolution of the system is ruled by 
the gravity due to the black hole and the scalar field 
only; that is, we assume our system occurs much after 
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the early universe stage, unlike previous studies related 
to primordial black holes growth through the accretion 
of a scalar field [HI; (ii) in most of cosmological models 
involving scalar fields, a first approximation of homogene- 
ity and isotropy is assumed 2, 4], however at local scales 
-for instance close to the event horizon of a black hole- 
wc consider such conditions do not hold, (iii) sometimes 
properties of pressure an density, and thus an equation 
of state, are attached to cosmological scalar fields |2|. [l8j . 
however in general at local scales scalar fields develop 
anisotropic pressures and time and space dependent re- 
lations between density and pressure (seejl9| for scalar 
field self-gravitating configurations and [20[ for the ac- 
cretion of a phantom scalar field onto a black hole) that 
we consider irrelevant in our analysis. 

This paper is organized as follows. In Sec. |TT] we de- 
scribe the system of equations ruling the evolution of the 
scalar field and the geometry, in Sec. IHII we describe the 
numerical methods used to solve the equations. In sec 
llVl we describe our results and in Sec. |V]we discuss our 
results and draw some conclusions. 
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Wirp ADAI = 0, 
j ij VjK - 8irf = 0, 



(4) 



where ^R is the scalar of curvature associated to 7^. In 
the equations Q|, the quantities padm, j l , Sij and 
S = 7 U Sij correspond to the local energy density, the 
momentum density, the spatial stress tensor and its trace 
respectively, measured by an Eulerian observer. These 
variables are obtained from the projection of the energy 
momentum tensor along the spacelike hypersurfaccs 
and along the normal direction to such hypcrsurfaces. 

We will now restrict to spherically symmetric black 
holes. In such case we consider the following ansatz for 
the 3- metric 7^ , the extrinsic curvature Kij and the shift 
vector (3 l : 



7ij 







j ee sin 2 9 



(5a) 



II. EVOLUTION EQUATIONS 

A. Evolution of the geometry 

In order to solve numerically the Einstein Field equa- 
tions G^ v = 87rT Mi , , where G^ u is the Einstein tensor and 
T^y is the energy momentum tensor, we use the 3+1 de- 
composition space-time and adopt the ADM system of 
evolution equations. We start up with the metric 

ds 2 = -(a 2 - Pi/3 l )dt 2 + 2/3 i dx i dt + - ,//.- '(/./•'. (1) 



Ki 






Kgg sin 2 9 



^ = (/3 r ,0,0) 



(5b) 



(5c) 



where the usual polar spherical topology with spatial co- 
ordinates x l — (r, 9, ip) is used. The evolution of the 
space-time geometry described in general by ^ and (O 
reduces in the present case to the following set of equa- 
tions 



where a is the lapse function, ft 1 are the components 
of the shift vector, 7^ are the components of the 3- 
metric of the hyper-surfaces that foliate the space-time 
and a; M = (t,x z ) are the coordinates of the space-time. 
All our calculations assume geometric units G = c = 1 . 

According to the ADM formulation of general relativ- 
ity, Einstein's equations split into evolution equations for 
the 3-metric 7^ and the extrinsic curvature Kij of the 
hypcrsurfaces 



dtjij = —2aKij + V,/3j + Vjft, 



dtKij = - ViVjCt + a (Rij + KK 



2K,,K. 



+ 4ira [(S - pADM)lij 



- 2 S^] 
KjtVi/3 1 , 



(2) 



(3) 



where Vi denotes the covariant derivative with respect to 
the 3-metric, Rij is the Ricci tensor related to the space- 
like hypcrsurfaces and K = 7 y iQj is the trace of the ex- 
trinsic curvature. In addition to the evolution equations, 
there arc the Hamiltonian and Momentum constraints 



dtjrr = —2aK rr + /3 r d r j rr + 2~f rr d r (3 r , 
daee = -2aK ee + /3 r d r jgg, 

(d r -f rr )(d r a) a( d r ^gg x 2 
2 



d t K„ 



-d rr a 
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2j rr jee 
/3 r d r K rr + 2K rr d r /3 r 



lee 



7 rr 

+ AlTa[(S - PADM)lrr - 25 rr 

c, {d r -ygg)(d r a) d rr Jgg 
o t Kee = a- 



2 r y rr 

(d r Jrr)(d r j6e) , 

a — "2 h a 



2 lr 



+ Aira[{S - PADAihee - 2S 



P T d r Kgg 

(6) 



These are the equations of the evolution of geometry 
that have to be evolved together with the matter field 
equations. 
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B. Evolution of the scalar field 

The matter held in the equations above corresponds to 
a scalar field described by the stress-energy tensor 

V = - \g,Jp' S *,s + m 2 3> 2 }, (7) 

whose evolution equation is given by the Bianchi identi- 
ties T^ u '.„ = 0, which reduces to the Klein-Gordon equa- 
tion 

□ $ - ?ti 2 $ = 0, (8) 

where the D'Alambertian operator for a general space- 
time is = ■j=d li [y/ II gg> 1 ''d v $]. Like in [IH, we re- 
strict to the mass-less case as a first approximation of 
the full study of the parameter space of the scalar field 
properties. In order to couple the evolution of the geom- 
etry and matter with the same evolution algorithm we 
write equation ([5]) for m = as a first order system of 
equations as suggested in plj : 

d t U = d r ( a y r ^) + aKU + (3 r d r U, 

d t $ = an + f3 r y, (9) 

where II = (<9 t $ — /3 r d r G>)/a and = <9 r $ are new 
first order variables. In terms of these new variables, 
the source terms in the ADM equations read: 

tepADM = \ (7 rr * 2 + n 2 ) , 
Anj r = -*II, 

47r5 = -i 7 rr * 2 + ^n 2 . (10) 

In summary, the evolution equations to be solved are (|6|) 
and ((9]) subject to the constraints ((4]). 

III. NUMERICAL METHODS 

We solve the system of equations as an initial value 
problem using a finite differences approximation on a sin- 
gle resolution uniform grid. A description of each stage 
of the solution is as follows. 



A. Initial Data 

In order to start up the evolution it is necessary to 
construct initial data consistent with Einstein's equa- 
tions, unlike the cases where the scalar filed is a test 



field [HI, [T5l - [T7| , where arbitrary initial scalar field pro- 
files are used on top of a given background. We simplify 
the constraints (j4|) by assuming the ansatz jgg — r 2 and 
K rr — _2M_i+m^j_ ^.jjjgjj re d uce the constraints flU) 

^ y/l+2M/r ' *-* 

to the following set of ordinary equations 

Q Irr A lrrKgg\ „7rr r , r , 
Orjrr = 1 — Irr n — ^ K rr Kgg 

r \ t J t 

+ 8irrj 2 r p ADM , 

dr Kee = r§^ + ^-47rr 2 lrr f. (11) 
Kee r 

We solve these equations once we have the explicit matter 
fields, for which we need to pro vide a scalar field profile 
similar to that described in [15| which corresponds to a 
spherical wave modulated by a Gaussian profile: 

$(t = , r)= A(- > 2 /^, 

r 

U(t = 0,r) = - — $(t = 0,r). (12) 
a 

With this information we calculate padm and j r at initial 
time using (|10[) . then substitute into equations (1111) and 
solve using a fourth order Runge-Kutta integrator. 

Unlike the problems where the space-time is assumed 
to remain fixed (e. g. [HI, [HI, the space-time is 
initially distorted with respect to an exact black hole ge- 
ometry. This fact makes difficult to filter ingoing and 
outgoing pure modes of the scalar field initially as de- 
scribed in [l5j . Therefore we do not attempt to have 
only ingoing or outgoing pure modes and will keep the 
presence of both. Since we want a mass of reference to 
compare with the growth of the black hole's horizon, we 
calculate the Misner-Sharp mass (see below) which con- 
tains the information -in terms of mass- of the ingoing 
and outgoing pulses. We will thus focus on the Misner- 
Sharp mass related to the ingoing pulse and define the 
ADM mass as the Misner-Sharp mass of the space-time 
minus the mass carried by the outgoing pulse. 

B. Evolution 

In order to solve the evolution equations (O and © we 
use the method of lines with a third order Runge-Kutta 
integrator and second order accurate spatial stencils. On 
the other hand, Bianchi identities guarantee that if the 
Hamiltonian and momentum constraints are satisfied at 
initial time, they would be satisfied further. We moni- 
tor that the constraints are truly being satisfied in the 
continuum limit using convergence tests. 

We avoid the singularity of the space-time at r = 
excising a chunk of the domain inside the event hori- 
zon 1221] . Thus we define the radial domain such that 
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r £ [nn,r ext } with < r in < r horizon and r ext as big as 
possible. Provided we use Eddington-Finkesltein slices 
to describe the black hole space-time, light cones at Tin 
are open and there is no need to impose any boundary 
conditions there. On the other hand, at r = jvw we ap- 
ply outgoing radiative boundary conditions [21| . Even 
though, in all our runs we make sure the exterior bound- 
ary is causally disconnected from the region where we 
measure physical quantities, in order to avoid the pol- 
lution from numerical errors reflected from the exterior 
boundary. 

We update the gauge functions a and (3 r using the 
conditions a/ '^-f r r + P r — 1 and ~2aKee + (3 r d r j96 = 
that guarantee that t + r is in ingoing null coordinate and 
that <9t76>0 = respectively plf . 

C. Diagnostics 

Since we want to study the accreted scalar field, it is 
useful to track the apparent horizon of the black hole, 
which provides an approximate growth rate of the black 
hole mass, which further indicates the correct accreted 
scalar field. We achieve this by locating the outermost 
marginally trapped surface, which obeys the condition 

= ^-^=0, (13) 

700 \Jlrr lee 

which is a procedure we can practice on the fly during the 
evolution Y2M. Then we define the radius at which this 
happens as the apparent horizon radius Tah = \flee~AH 
in terms of the areal radius. Then we estimate the ap- 
parent horizon mass using Mah — tah/^- Since our 
results converge with second order due to the discretiza- 
tion along the spatial direction, the final apparent hori- 
zon mass is a Richardson extrapolation of our results 
using the two finest resolutions in all our calculations. 

We also locate the event horizon because it is a gauge 
independent 3-surface unlike the apparent horizon which 
is a gauge dependent 2-sruface. The reason to locate the 
event horizon is that we want to make sure the event 
horizon lies outside the apparent horizon and close to 
it after the accretion of the scalar field has finished. We 
locate the event horizon as the attractor surface of inward 
null rays triggered from the future and from outside the 
black hole. 

In order to make sure that the horizon growth is consis- 
tent with the mass of the space-time including the contri- 
bution of the scalar field, we measure the Misner-Sharp 
mass [23| defined for our space-time as 

M '^(l-I<^ + ^) (14, 
2 V 4 IrrlOO 700 J 

which allows us to estimate the ADM mass Madm = 
limj._j.oo Mms- In this way we can diagnose whether or 



not the final mass of the black hole's horizon corresponds 
to Madm- 

Finally we also track the constraints and practice 
the required convergence tests they must satisfy. All our 
results assume the radial and time coordinate are written 
in units of M, where M is the initial mass of the black 
hole's event horizon. 



D. Example 

Now we present in detail one of our simulations. In Fig. 
[1] we illustrate how we estimate the mass of the space- 
time. Some snapshots are shown for the scalar field and 
for the Misner-Sharp mass. We show on the left panels 
how the initial pulse splits into two pulses, one moving 
toward the black hole and the other moving outwards. 
On the right panels wc show the Misner-Sharp mass and 
illustrate how the outgoing pulse carries a given amount 
of the space-time mass. The mass measured once the sys- 
tem has stabilized is what we call the ADM mass Madm, 
around t ~ 75. This would be the mass of a space-time 
if it would only containing the ingoing pulse. Finally in 
Fig. [I] we present the convergence test of the Li norm of 
the Hamiltonian constraint in order to illustrate how all 
our runs behave. 

In Fig. [5] we show other of our simulations showing 
both the apparent and the event horizons. What we want 
to illustrate is that on the one hand the apparent horizon 
lies always inside the event horizon, which is consistent 
with the energy conditions of the space-time and also, 
that once the black hole has accreted the incoming scalar 
field, both horizons coincide, which authorizes us to use 
the apparent horizon as the surface to monitor the mass 
of the black hole in time. 



IV. RESULTS 

We split the parameter space in terms of the wave num- 
ber k, which is one of the properties of the propagation 
speed of the wave packet. In fact it was found for a fixed 
background that for k < 1 the absorption rates dimin- 
ish. Thus we choose four values of k = 0, 0.5, 1, 2, for 
which wc choose a rather spread set of parameter values 
A = 0.5, 0.8 and a = 1, 5, 10. With these parameters 
we sample different scalar field contributions to the en- 
ergy density of the space-time, length scales and number 
of nodes of the wave function. In [l5j the parameter a is 
related to the variance on k provided the wave number 
lies around ko such that ((fc — ko) 2 ) = 1/cr 2 and full ab- 
sorption was found for k > 1, and only partial absorption 
in other cases. Wc choose two values of k smaller than 
the threshold k = 1 and two bi gge r in order to confirm 
or contradict the predictions in [151 ] . 

Our results are shown in Figs. [31 HI E] and [HI that 
correspond to the wave numbers k = 0, k = 0.5, k = 
1 and k ~ 2 respectively, for the different values of A 
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k=2, A=0.8, o=10 




40 
f 



FIG. 1. We show one of the cases of our sample in detail, the 
one corresponding to k = 2, A = 0.8 and a — 10. We show 
snapshots of the pulse moving toward the black hole. As men- 
tioned in the text, the initial pulse splits into two, one moving 
inward and the other outward. We show on the right sides the 
Misner-Sharp mass for the same snapshots and show that it 
approaches a constant value that includes only the incoming 
pulse contribution. The Misner-Sharp mass becomes a con- 
stant after the two pulses carry their respective mass-energy 
contributions. This constant value is what we call our Mad m 
and is free of the contribution of the outgoing pulse. Finally 
we present as a sample of all our convergence tests the L2 
norm of the Hamiltonian constraint. The resolutions used in 
our runs are Ari = 0.0125M (continuous line), Ar2 = An/2 
(dashed line) and Ar$ = Ar?/2 (dotted line). 



Apparent Horizon 




Event Horizon 



1.5 2 2.5 3 3.5 4 
r 

FIG. 2. We show one of the cases of our sample in detail, the 
one corresponding to k = 1, A = 0.8 and cr = 5. This plot 
shows how the apparent horizon and event horizon coincide 
both at initial time and after the black hole has stabilized 
after the accretion of the incoming scalar field pulse. We also 
show a set of outgoing null rays, or equivalently ingoing past 
directed null rays, that converge to a surface that happens to 
be the event horizon. The event horizon is only tracked up to 
t ~ 80, after which our null rays diverge. 



and cr. In these figures we show the apparent horizon 
and Madm- These parameters are organized in such a 
way that we explore the suggestions in related to 
the absorption of the scalar field in terms of the wave 
number, amplitude and width of the wave packet. In all 
our runs we use rg = 40, which is far enough from the 
horizon. 

Our first results confirm that for k > 1 full absorption 
is observed independently of the value of a (Figs. [5j[6]). 

In the fixed space-time background case it was found 
that for values of k smaller than 1 and big values of a not 
all the incoming scalar field was absorbed 15j. However 
as shown in Figs. [3J H] when k = 0.5 full absorption was 
found in all cases and the only case where we find non- 
full (partial) absorption is that of k = 0. Specifically for 
a < 2, that is, smaller than the Schwarzschild radius, 
full absorption is found, whereas for a = 5 for the two 
values of A we use, only ~ 91% of the ingoing scalar field 
is absorbed and for a = 10 we found that ~ 65% of the 
ingoing scalar field is absorbed. For comparison, in the 
fixed background case, it was found that the absorption 
could be as small as 50% Hal. 



V. DISCUSSION AND CONCLUSIONS 

We have tracked the spherically symmetric non-linear 
evolution of a massless scalar field being accreted into 
a Schwarzschild black hole. The aim of the work is to 
determine whether or not a scalar field is reflected or 
absorbed in terms of its wave number or wave packet 
width. As a first approximation we have considered the 
field to be massless. 
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k=0, A=0.5, o=l k=0, A=0.8, o=l >W0.5, o=l /(=0.5, o=l 




FIG. 3. We show the time evolution of the apparent horizon 
mass and the Madm for the case k = 0. The fact that the 
apparent horizon mass (continuous line) grows up to the value 
Mad m (dashed line) indicates that the whole incoming scalar 
field energy density has been absorbed, which happens only 
for (7=1, smaller than the Schwarzschild radius rs = 2, 
whereas for the bigger values a = 5, 10 a fraction of the 
scalar field is not being absorbed. The proportion of the field 
depends only on the value of a, and we found that for a = 5 
92% of the energy density of the scalar field has been absorbed 
and for a = 10 67%. 



With our results we generalize the analysis of the reflec- 
tion and absorption of a massless scalar field in order to 
bound the possibility that black holes, like supermassive 
black holes could potentially absorb scalar fields that are 
being used as ingredients of nowadays cosmological mod- 
els like scalar field dark matter and quintessence. Our 
generalization consists in adding the evolution of the ge- 
ometry of space-time, unlike previous analyses where it 
remains fixed. 

Our results confirm some of the previous predictions 
related to fixed background analyzes, specifically, we con- 
firm that for spherical wave profiles of the scalar field with 
k = not all the incoming energy density is absorbed. 
Instead, we found the tendency to have less absorption 
when the initial wave packet width is bigger. 

For all the other values of k used, we found full absorp- 
tion. Nevertheless, a detailed study would include other 
values of k closer to 0. 

In the context of the scalar field dark matter and dark 
energy models consisting of scalar fields, our results indi- 
cate that the evolution of the geometry allows the absorp- 
tion of the total amount of the scalar field that is radially 
directed toward the black hole within a short time scale 



FIG. 4. We show the time evolution of the apparent horizon 
mass (continuous line) and the Madm (dashed line) for the 
case k = 0.5. Unlike the case k — 0, in the present case full 
absorption was found in all cases. 



k=1 , A=0.5. a=\ k-1, A-0.8, obI 




( f 



FIG. 5. We show the time evolution of the apparent horizon 
mass (continuous line) and the Madm (dashed line) for the 
case k — 1. Also full absorption was found in all cases. 
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k-S, A=0.5. o=l k.2, A-0.8, o=l 



( I 
k=2, /3=0.5, <J=5 k.2, A=0.8, o=5 




FIG. 6. We show the time evolution of the apparent horizon 
mass (continuous line) and the Madm (dashed line) for the 
case k = 2. Also full absorption was found in all cases. 
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for large wave numbers, whereas, for k ~ a fraction of 
the scalar field remains outside the black hole. 

At this point we have only performed our analyses us- 
ing pretty general initial scalar field profiles. Nevertheless 
it would be interesting to study evolutions with quasi- 
stationary initial data (on a fixed background) like those 
proposed in [lj} and study their evolution on an evolving 
space-time geometry. 

Also, if specific models of dark matter or dark energy 
are to be studied, it is necessary to include the mass 
term of the scalar field and less restrictive symmetries, 
like rotation of the scalar field near the black hole and 
Kerr black holes. 
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